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ABSTRACT 
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An investigation is made of disk accretion of matter onto a rotating star with an 
! ahgned dipole magnetic field. A new aspect of this work is that when the angular velocity 
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of the star and disk differ substantially we argue that the B field linking the star and disk 
rapidly inflates to give regions of open field lines extending from the polar caps of the 
^ ! star and from the disk. The open field line region of the disk leads to the possibility of 
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magnetically driven outflows. An analysis is made of the outflows and their back affect on 
the disk structure assuming an "ct" turbulent viscosity model for the disk and a magnetic 
^ . diffusivity comparable to this viscosity. The outflows are found to extend over a range of 



radial distances inward to a distance close to rto, which is the distance of the maximum 



O ' of the angular rotation rate of the disk. We find that rto depends on the star's magnetic 
c/3 . moment, the accretion rate, and the disk's magnetic diffusivity. The outflow regime is 

a ■ 

accompanied in general by a spin-up of the rotation rate of the star. When rto exceeds the 
^ ' star's corotation radius = (GM/cu^)^, we argue that outflow solutions do not occur, 
but instead that "magnetic braking" of the star by the disk due to field-line twisting occurs 
in the vicinity of rcr- The magnetic braking solutions can give spin- up or spin-down (or 
no spin change) of the star depending mainly on the star's magnetic moment and the 
mass accretion rate. For a system with rto comparable to r^, bimodal behavior is possible 
where extraneous perturbations (for example, intermittency of a, B field fiux introduced 
from the companion star, or variations in the mass accretion rate) cause the system to fiip 
between spin-up (with outfiows, rto < fcr) and spin-down (or spin-up) (with no outfiows, 
rto > rcr)- 

* On leave from Department of Applied Physics, Cornell University, Ithaca, NY 14853. 
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I. INTRODUCTION 

The problem of matter accretion on to magnetized stars has been of continued interest 
since the discovery of X-ray pulsars in binary systems (Giacconi et al. 1971; Schreier et 
al. 1972; Tananbaum et al. 1972). The X-ray pulsars were interpreted as magnetized 
neutron stars, accreting matter from a companion star (Pringle and Rees 1972; Davidson 
and Ostriker 1973; Lamb, Pethick, and Pines 1973; Rappaport and Joss 1977). Of the 
more than 40 of known X-ray pulsars, most are in high-mass {Mfot > 15) binary systems, 
while only few are in low-mass {Mfot < 3) binary systems (see, for example, review by 
Nagase 1989). 

From the early observations, many X-ray pulsars were found to show secular decreases 
of pulse period (spin-up of the neutron star rotation) with relatively large rates of change 
of —P/P = 10~^ — 10~^ yr"^ (Gursky and Schreier 1975; Schreier and Fabbiano 1976; 
Rappaport and Joss 1983). The long-term monitoring of X-ray pulsars during the last two 
decades have shown a wide variety of pulse-period evolution on different time-scales. Some 
X-ray pulsars show a steady spin-up evolution during many years which rapidly changes 
to a steady spin-down evolution (e.g., GX SMC X-1, 4U 1626-67). Others show the 
wavy variations of P on time-scales of a few years on a background of systematic spin-up 
or spin-down (e.g.. Gen X-3, Her X-1, Vela X-1) (Nagase 1989; Sheffer et al. 1992; Bildsten 
1993). 

Short-term fluctuations of pulse-period on time-scales of days to months, including 
clear evidence of spin-down episodes, were found for a number of sources, for example. Her 
X-1 (Giacconi 1974), Gen X-3 (Fabbiano and Schreier 1977), Vela X-1 (Nagase et al. 1984). 
Detailed observations of Vela X-1 with Hakucho and Tenma satellites, and with HEAO-1 
have revealed fluctuations of large amplitude with time-scales about 3-10 days, which were 
estimated to be random in both amplitude and sign (Nagase 1981). BATSE observations 
have shown that in some cases (Gen X-3, Her X-1) the timescale of fluctuations is less than 
one day with the change of sign from spin-up to spin-down and vice- versa (Bildsten 1993). 
Note that some X-ray pulsars (usually the long-period ones) show no evidence of regular 
spin- up or spin-down (for example, 4U 0115-1-63, GX 301-2). 

In some cases, like Her X-1, the existence of an accretion disk is indicated from various 
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observations in the X-ray and optical bands (for example, Middleditch and Nelson 1976; 
Bisnovatyi-Kogan et al. 1977; Middleditch, Puetter, and Pennypacker 1985). In several 
other cases, such as Cen X-3, SMC X-1, GX 1+4, and 4U 1626-67, there are also different 
kinds of evidence, that matter from companion star forms a disk (for example, Tjemkes 
et al. 1986; Nagase 1989). Most of these are short-period pulsars with period about 
several seconds or less (excluding GX 1+4 with period 122 s.). However, most of long- 
period pulsars, are supposed to be powered by the capture of stellar wind from companion. 
However if the companion is the star of Be type with slow rate of outflow, then it is probable 
that once again the disk first forms, opposite to the case when the companion star is of 
0-B kind with high speed of matter outflow (see review of Nagase 1989). Illustrative data 
for two pulsars are shown in Figure 1. 

Some aspects of matter accretion by magnetized neutron stars were considered for the 
first time by Bisnovatyi-Kogan and Friedman (1970), and by Shvartsman (1971). Ideas on 
the spinning-up of pulsars by accretion were first proposed in the case of disk-fed pulsars by 
Pringle and Rees (1972), Lamb, Pethick and Pines (1973), and Lynden-Bell and Pringle 
(1974), and for wind-fed pulsars by Davidson and Ostriker (1973), and lUarionov and 
Sunyaev (1975) (see also Lipunov 1993). The spin-down of a rapidly rotating magnetized 
star due to "propellar" action of the magnctosphcrc was suggested by Shvartsman (1971) 
and lUarinov and Sunyaev (1975), but a detailed model has not been worked out. Pringle 
and Rees (1972) considered matter accretion from a viscous disk to a magnetized rotating 
neutron star. Their theory considered the idealized case of an aligned rotator where the 
magnetic moment of pulsar is aligned parallel or anti-parallel with the pulsar's spin-axis 
and the normal to the plane of the accretion disk. They supposed that accreting matter 
stops at the point where the magnetic pressure of the magnetosphere becomes equal to 
pressure of matter. At this distance, matter transfers its angular momentum to the star 
through a thin boundary layer. Another idea regarding the width of the transition zone 
between the unperturbed accretion disk and magnetosphere, was considered in scries of 
papers by Ghosh and Lamb (1978, 1979a,b) (hereafter denoted GL). They assumed that 
magnetic field of the neutron star threads the accretion disk at different radii in a broad 
"transition" zone, in spite of predominance of the disk stress compared with magnetic 
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stress in this region. A schematic diagram of their model is shown in Figure 2. Their work 
is based on the idea of anomalous resistivity connected with the supposed fast reconnection 
of toroidal component of the magnetic field lines across the disk. 

Here, we propose an essentially different picture for disk accretion onto a star with 
an aligned dipole magnetic field. The idea is based on the fact that when there is a large 
difference between the angular velocity of the star and that of the disk, the magnetic field 
lines threading the star and the disk undergo a rapid inflation so that the fleld becomes open 
with separate regions of field lines extending outward from both the star and the disk. As a 
result the magnetosphere consists of an open field line region far from the star and a closed 
region approximately corotating close to the star. This is shown schematically in Figure 
3. Our model does not assume an anomalous resistivity, but rather a turbulent magnetic 
diffusivity of the disk comparable to the turbulent a viscosity proposed by Shakura (1973). 
This corresponds to a resistivity about 10^ times smaller than that assumed by Ghosh and 
Lamb (1978,1979a,b). Campbell (1992) has earlier discussed closed field line models of 
disk accretion onto an aligned dipole star assuming a magnetic diffusivity comparable to 
the turbulent a viscosity. 

Here, we study the matter flow in the disk taking into account the open magnetic 
field line region of the magnetosphere. The existence of an open magnetic field region of 
the disk leads to the possibility of magnetically driven outfiows. We analyze the magne- 
tohydrodynamic (MHD) outflows and their back influence on the disk using the work on 
MHD outflows and magnetized disks by Lovelace, Berk, and Contopoulos (1991) (denoted 
LBC), Lovelace, Romanova, and Contopoulos (1993) (LRC), and Lovelace, Romanova, and 
Newman (1994) (LRN). 

In §n we discuss the basic equations. In §IIA we argue that part of the star/disk 
field configuration is open. In §IIB we give the basic equations for a magnetized viscous 
accretion disk. In §IIC we discuss magnetically driven outflows and in §IID, necessary 
conditon for these outflows. In §IIE, we give the main results for the case of outflows 
where in general star spins up. In §IIF, we discuss the magnetic braking of the star by 
the disk due to fleld line twisting, which occurs when there are no outflows, and show that 
the star can either spin-down or spin up depending on the mass accretion rate and the 
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magnetic moment of the star. In §111 we discuss the conclusions of this work. 
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II. THEORY 

The basic equations for an assumed stationary configuration of plasma are: 



V • (pv) = , V X B 




V-B = 0, VxE = 0, 



J = c7e(E + V X B/c) , p(v • V)v = -Vp + pg + -J X B + 



(1) 



Here, v is the flow velocity, p is the density, cTg is the effective electrical conductivity, F^"^ is 
the viscous force density, p = pkBT/m is the gas pressure (with ks the Boltzmann constant, 
m the mean particle mass, and with the radiation pressure assumed negligible), and g is 
the gravitational acceleration. Outside of the disk, dissipative effects are considered to be 
negligible (cTe ^ oo, F'^'® = 0, etc.). We neglect the self-gravity of the disk and relativistic 
effects so that g = -V$g with $g = -GM/(r^ + 2^)^/^, where M is the mass of the 
central star. Equations (1) are supplemented later by a equation for the conservation of 
energy in the disk. 

A general axisymmetric B- field can be written as B = Bp + 0-B<^, where Bp = V x 
{(j)^/r) is the poloidal field, is the toroidal field, and ^ is the flux function (see, 
for example, Mestel 1968, or Lovelace et al. 1986). We use a non-rotating cylindrical 
coordinate system so that Bp = (S^, 0, Bz). Notice that \l/(r, z) = const, labels a poloidal 
held line, (Bp • V)^' = 0, or a flux-surface if the poloidal held line is rotated about the 
z-axis. For the present problem, the Bp fleld can be represented as = + with 
the star's fleld-assumed to be a dipole = p.r'^{r'^ + z^)~i and with ^' due to the 
non-stellar toroidal currents. 

A. Inflation and Opening of Coronal B-Field 

Figure 4 shows a sketch of the poloidal projections of two nearby fleld lines connecting 
the star and the disk. Because d/dt = 0, the E fleld is electrostatic and the poloidal plane 
line integral of E around the loop, 1— >2— >3^4^1, is zero. Because of axisymmetry 
and the fact that E + vxB/c = outside of the disk, the line integrals of E along the 
two curved segments in Figure 2 vanish separately. That is, the electrostatic potential 
is a constant on any given flux surface (see, for example, Lovelace et al. 1986). The 
potential difference between the points (1,2) on the star's surface is — 5ri2 • E* = oj^^S'^/c, 
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with 6^ = — where we have assumed that the star is perfectly conducting [E* = 
— (v X B)*/c], and where a;* is the angular rotation rate of the star. Thus, we have 
drs4,{Er)d = — c<;*5\l//c = —uj^rBz{r, 0)/c, where the path element 5r34 is in the mid plane 
of the disk, and the "d" subscript indicates that the quantity is evaluated in the disk. In 
turn, we have Er\d = {Jr/ <^e — Vci)Bz/ c)d- Prom Ampere's law, Jr = —{c/4:7r){dB(j)/dz), the 
fact that B(i) is necessarily an odd function of z, and the approximation dB^/dz = {B^)h/h, 
we have Jr ~ —{c/An){B(f,)h/h, where h is the half-thickness of the disk, and the notation 
(■ • ■)h indicates that the quantity is evaluated aX z — h. Combining these results gives for 
stationary conditions 

hv 

{B^)h{r) = - — [uJd{r)-u.]B,{r) , (2) 

where rjt = (? I (47rcre) is the magnetic diffusivity of the disk, and uj^ is the angular rotation 
rate of the disk. The fractional variation in B^ between the midplane and surface of the 
disk is negligible [0(/i/r)] if the disk is thin in the sense that hjr <^ 1. For the conditions 
of interest here, we show in §IIB that h/r < Cg/vx, where Cg = {kBT/m)^ the Newtonian 
sound speed based on the mid- plane temperature of the disk, and vk = (GM/r)^ is the 
Keplerian velocity at r. Campbell (1992) has independently derived equation (2) but does 
not place any limit on the ratio \B(j)/Bz\ (see discussion following equation 3). 

In contrast with the work of Ghosh and Lamb (1978, 1979a,b), we assume that the 
magnetic diffusivity of the disk rjt is of the order of magnitude of the disk's effective viscosity 
(see Bisnovatyi-Kogan and Ruzmaikin 1976; Parker 1979; Campbell 1992). Additionally, 
we assume that the disk viscosity is the turbulent viscosity as formulated by Shakura (1973) 
and Shakura and Sunyaev (1973), Ut = {2/3)aCsh, where a is a dimensionless quantity less 
than unity. In this work, we assume that a is in the range 10~^ to 10~^. The possible 
contribution to the turbulent momentum flux due to small-scale magnetic field fluctuations 
is assumed included in a (Eardley and Lightman 1975; Coroniti 1981; Balbus and Hawley 
1992; Kaisig, Tajima, and Lovelace 1992). We write % = Dut where D is dimensionless 
and of order unity. 

With rjt the turbulent diffusivity, equation (2) gives \{B(j,)h/Bz\ — {3/2)r\u>d{r) — 
(jJ:^\/{aDcs). This ratio is a measure of the twist of the B held between the star and the 
disk. Por the outer part of the disk a;d(r) is expected to be close to the Keplerian value 
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(jJk = (GM/r^)^ . Thus, it would appear that the twist |(S<^)/i/5z| can be very much 
larger than unity, say, > 10. We argue here that such large values of the twist do not 
occur. 

For a discussion of the limitation on the twist \{Bfj))h/ Bg\ we consider that the plasma 
outside of the disk is force-free, J x B = 0, which is the coronal plasma limit of Gold 
and Hoyle (1960). This limit is valid under conditions where the kinetic energy density of 
the plasma is much less than B^/Stt. The general response of force-free coronal magnetic 
field loops to stress (differential twisting) applied to the loop foot points (in the solar 
photosphere) has been studied intensely by Aly (1984, 1991), Sturrock (1991), and Porter, 
Klimchuk, and Sturrock (1992). The conclusion of these studies is that a closed field loop 
with a small twist evolves into an open field-line configuration as the twist is increased. 
The related problem of the twisting of a force free magnetic field configuration with foot 
points at different radii in a Keplerian disk has been studied by Newman, Newman and 
Lovelace (1992), and Lynden-Bell and Boily (1994). The twisting due to the differential 
rotation of the disk acts to increase the total magnetic energy of the coronal field and this 
in turn acts to "inflate" the field as shown in Figure 5. 

An alternative measure of the field twist between the star and the disk is simply the 
difference in the azimuthal location, A(p, of the stellar and disk foot points of a given 
fiux tube. We have rdcp/dlp — Brj,{r, z) /\'Bp{r, z)\, where dip is the length element along 
a poloidal field line ^^(r, ^) = const. Because rB^{r.,z) depends only on in a force-free 
plasma, we have 

/disk jj 
tar ^ ■ 

For small A(p (say, < 1), the poloidal field is close to that of a dipole, and the integral 
gives A(j) = Ct[{B^)h/ Bz]d with Ct = 8/15. For increasing A^, the value of Ct increases 
because of the longer path length and the smaller value of iV^'j. The present situation is 
analogous to a current-carrying plasma column with a longitudinal B field. For the plasma 
column there is the well-known Kruskal-Shafranov stability condition (see, for example, 
Bateman 1980) on the twist Acp of the field (at the column's surface) over the length of 
the column. This condition, Acj) < 2tt, is required for stability against symmetry changing 
kink perturbations. A related, non-axisymmetric instability may accompany the above 
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mentioned inflation and opening of the star/disk field in response to increasing A0. 

We propose that there is a definite upper hmit on the twist A(/) of any field lines 
connecting the star and the disk. This limit implies a corresponding limit on \{B(f))h/ Bz\d 
in equation (2). If is larger than this limit, then the field is assumed to be open as 
indicated in Figure (3). Note that when the field threading the disk is open equation (2) 
does not apply. 

B. Basic Equations for the Disk 

The flow velocity in the disk is v(r) = — •u(r)f + t'<^(r)(/), where u is the accretion speed. 

h 

The surface density is cr = J dz p{r, z). The main disk equations are: 

-h 

1. Mass conservation: 

M = 27rrau, — const. (4) 

We assume that the change of M with r due to possible MHD outflows is small. 

2. Radial force balance: 

-h 

3. Angular momentum conservation: 

F - 2 , 'T^^t du 

F = r 0J+ — , (6) 

u dr 

where a; = v<^/r. The term —r^{B(j)Bg)h represents the outflow of angular momentum 
from the ±z surfaces of the disk (that is, a torque on the disk). We assume that the 
angular momentum carried by the matter of the outflow is small compared with that 
carried by the held. If the B held at r is open as discussed in §IIA, then the angular 
momentum outflow from the disk is carried to infinity. On the other hand, if the field 
is closed, then the angular momentum outfiow from (or infiow to) the disk is carried 
by the coronal B field to (or from) the star. 

4. Energy conservation: 



-h 



The first term is the viscous dissipation, the second is the Ohmic dissipation, while 
the third term is the power per unit area carried off by radiation (in the ±z directions) 
from the disk which is assumed optically thick. Here, k is the opacity assumed due 
to electron scattering, a and as = ac/4 are the usual radiation constants, T is the 
midplane temperature of the disk, and Tg// is disk's effective surface temperature. 

5. Conservation of magnetic flux for a general time-dependent disk: 

where the generally small radial diffusion of the B^g field has been neglected (see LRN). 
The first term inside the square brackets represents the advection of the field while 
the second term represents the diffusive drift. In a stationary state, 

Pr{r) = = — . (86) 

Bz Vt 

6. Vertical hydrostatic equilibrium: 

The 2;-component of the Navier-Stokes equation (1) gives the condition for vertical 
hydrostatic balance which can be written as 

''^%(.(^)-(^)^ = 0, (9) 

\rj vk 

where is the Newtonian sound speed based on the midplane temperature of the 
disk, and b = r (-B^)^ + (-B^)^ / (^Anav]^) (Wang, Lovelace, and Sulkanen 1990). 
Radiation pressure is negligible for the conditions of interest. For b <^ 2cs/vk, this 
equation gives the well-known relation h/r = Cs/vk (Shakura and Sunyaev 1973), 
while for b ^ 2cs/vk it gives h/r = b~^(cs/vKY which is smaller than Cg/vK owing to 
the compressive effect of the magnetic field external to the disk (Wang et al. 1990) It is 
useful to write b = ep^, where e = rSf (r, 0)/(27Tavj^) and = [(Br)l+iB^)l]/(2B^). 
In §IID, we discuss that /3 = 0{1) in order to have outflows from the disk, and 
/? < 1 with no outflows. Thus, the Alfven speed in the midplane of the disk is 
va = vk^^ {h/r)^ . We may term the magnetic held as weak if e < Cs/vk ~ h/r. In 
this limit, va/cs ^ (er/h)^ < 1, and the magnetic compression of the disk is always 
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small. In the opposite, strong field limit, e > Cg/vx- If at the same time (3 = 0(1), 
then the disk is magnetically compressed, and for e ^ Cs/vk, va/cs = 

The following subsections involve applications of equations (4) - (8). Note that with 
the magnetic field terms neglected these equations give the Shakura-Sunyaev (1973) solu- 
tion for region "6" . 

C. Magnetically Driven Outflovirs 

Consider the outer region of the disk where the value of the field twist A(p given 
by equation (3) [or \{B(i))h/Bg\ of equation (2)] is larger than a critical value. In this 
region the B field threading the disk is open as discussed in §IIA. The angular momentum 
flux carried by the disk is MF. If the angular rotation of the disk is approximately 
Kcplerian, u{r) ^ UK{r), then the viscous transport contribution to F is r'^ {u / u) {duj / dr) ~ 
— {3/2)r'^uj[h/ {(3rr)]. As discussed in detail below a necessary condition for MHD outflows 
is that Pr = Consequently, the viscous contribution to F in equation (6) is smaller 

than the bulk transport term by the small factor h/r « 1, and F cur^. 

In equation (6) we therefore have d{MF)/dr ^ Mujkt/I. In equation (6) we let 



{B^)h = -(3^B^{r,0). Studies of MHD outflows (Blandford and Payne 1982, LBC, LRC) 
indicate that = const. ^ 0{1). As a result, equation (6) implies that 



where the w subscript indicates the wind region of the disk. Of course at a sufficiently 
small radius, denoted r^i , the outfiow ceases, and the dependence of Bg{r,0) reverts to 
approximately the stellar dipole field, ///r^. As a consistency condition we must have 



The transition between the dipole and outflow fleld dependences is handled by letting 
Bz{r,0) = {pi,/r^)g{r) + {k/r^/^)\l — g{r)], where, for example, g{r) = {1 + exp[(r — 
ri„j)/Ar]}"^ and /^r/ry,i -C 1. 

Equations (4) - (9) can be used to obtain the disk parameters in the region of outflow. 
As shown below the viscous dissipation is much larger than the Ohmic, and equation (7) 



[i?,(r,0)]^ = k/r'''\ with k = [M{GM)-^ / {2(3^)]^ , 



(10) 




(11) 
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then gives T = [81M^{GMy/^K/ {1287r^aDp^ac)]^r-^ . For illustrative values, M = 
IMq, M = 10^ Vs ~ 1-6 X lO-^Mg/yr, a = 0.1, pr = 0.58, and i:> = 1, we find at 
the representative distance r = lO^cm that T ^ 0.44 x 10^ K, vk = 1.15 x lO^cm/s, 
Cs = 7.6 X 10*^cm/s, u = 0.88 x 10*^cm/s, a = 180g/cm2, and h/r = 6.6 x 10"^. The 
ratio of the ohmic dissipation to that due to viscosity is: (4/27)q;D/3^ -C 1. The ratio 
of the viscous dissipation to the power output of the outflows (the ± z Poynting fluxes) 
(per unit area of the disk) is: {9/2){h/r){PrD)~^ <S 1. Thus, most of the accretion power 
for r > r^i goes into the outflows. However, the fraction of the total accretion power in 
outflows is small if r^i ^ r*, where r* is the star's radius (~ 10^ cm for a neutron star). 
Note that M = IQ-'^^g/s and M = IMq correspond to a total accretion power or luminosity 
Lq = GMM/r^ fa 1.33 x 1037(erg/s)(106cm/r*). 

Conservation of the poloidal flux threading the disk implies that there is an outer 
radius, denoted r^o^ of the outflow from the disk. That is, we have f^""" rdr(k/r^'^) = 
r'^° rdriu/r^). Thus, the outer radius is r^o ~ fwiif^)~'^^^i where failforl — £<Cl 
and / = 3/4 for f <^ 1. For r > Vwoi Bz{r, 0) rapidly approaches zero. 

Equation (9) gives the disk thickness. From equation (10), we have e = rB'^/{27rav%) 
= {u/cs){cs/vK)/{2(3(f,). The magnetic compression of the disk is small if e < Cs/vk or 
equivalently if u/cg < 1. Our numerical solutions have u/cs <1 in the region of outflow. 
D. Necessary Condition for MHD Outflows 

A necessary condition on for MHD outflows from the disk can be obtained by 
considering the net force J^p on a fluid particle in the direction of its poloidal motion 
above the disk. Notice that for stationary flows, the poloidal flow velocity Vp is parallel 
to Bp owing to the assumed axisymmetry and perfect conductivity. However, in general 
V X B so that the fluid particles do not move like "beads on a wire" . Above, but close 
to the disk {h < z <^ r),we assume that the poloidal fleld lines are approximately straight. 
Thus, the poloidal position of a fluid particle above the disk is r = {vo + Ssin{6))r + 
Scos{6)z, where 5" is the distance along the path from the starting value So = h/cos{9) 
at a radius Tq, tan{6) = {Br)h/Bz = P^ and z = Scos{9). The effective potential is 
U{S) = -(l/2)a;2(ro + Ssin{9))^ - GM/\r\ (LRN), and the force in the direction of the 
particle's poloidal motion is J^p = —dU /dS. For distances z much less than the distance 
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to the Alfven point, Uq ~ 0JK{fo)- In this way we find 




(12) 



for /i ~ z <S r, where we assume > 0. 

The slow magnetosonic point of the outflow occurs where — Q (LRC) at the distance 



The factor 1 — [uj/uk)"^ can be obtained from the radial force balance equation (5). The 
radial pressure force is small compared with the magnetic force for the conditions we 
consider, q;/?^/(3/?<^) > h/r < 1. Thus, 



For Tyji < r < r-ujoi the quantity in brackets is simply q:/(3/3<^) owing to equation (10). 
Note that the minimum of /?^/(3/?^ — 1) is 1/2 at /J^ = 1. If the gas near the surface of 
the disk {h ^ z < Zg) is assumed isothermal with temperature Ti <^ T and sound speed 
Csi = (/csTi/m)^, then the density at the slow magnetosonic point ps can be obtained 
using the MHD form of Bernoulli's equation (LRC). At the slow magnetosonic point, 
the poloidal flow speed is Vp = Csi|Bp|/|B| = Vsm, the slow magnetosonic speed, and 
Ps = Phexp{— 1/2 — [U{zs) — t/(^)]/Csi}, where U is now regarded as a function of z. In 
turn, the mass flux density from the +z surface of the disk is psVsmCos{0) = psCsi-Bz/|B|. 

The consistency of the magnetically driven outflow solutions requires Zg ~ h. If this 
were not the case, the outflows would be matter rather than field dominated near the 
disk. There are two possible regimes having Zs ~ h: one is with aD/{?tj3(j>) < 1 [or 
aD/(3Pcf>) < 1] and Pr ~ 1/^3; another is with aD/{<oj3^) > 1 and /^^ > 1. In the 
present work, we consider the first regime which is consistent in the following respect: 
The magnetic pinching force on the inner part of the outflow (r < Two) increases as the 
mass flux in the outflow increases (LBC). In turn, an increase in the pinching force acts 
to decrease Pr = {Br)h/Bg thereby increasing Zg/h and decreasing the mass flux in the 
outflow. 




(13) 




(14) 
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In addition to the condition Zg ~ the slow magnetosonic point must be at a distance 
Zs not much larger than /i, say 2/i, in order for the outflow to have a non-zero mass flux. At 
the inner radius of the outflow, r^i , the S^(r, 0) field in the disk is /x/ (2r^j) as discussed 
in §IIC. We then have 



aDjj? 



24M(G'M)2 



0.75 X lO^cm 



VO.l/ VlOSOGcm^/ V M )\m) 



(15) 

from equation (14). Note that from equation (11), we get S ~ [q;D/(48/3<^)]^ or £ ~ 0.046 
and Tyjo ~ 90r^yj for a = 0.1, D — 1 and ~ 1. 
E. MHD Outflows and Spin-Up {uo < Vcr) 

We have numerically integrated equations (4) - (9) assuming magnetically driven 
outflows from the disk as discussed in §§IIC & D. We integrate the equations inward 
starting from a large radial distance r < ryjo ■ The inner radius of the region of outflow is 
determined as indicated in §IID. For r < Vwi , dF/dr = 0, and the solutions in this range 
of r all exhibit a "turn-over radius", rto j where dw/dr = 0. This "turn over" results from 
the radially outward magnetic force in equation (5) which becomes stronger as r decreases. 
In the region close to the turnover, the magnetic compression of the disk becomes strong 
in that e = 0(1). Figure 6 shows a representative solution. From a least squares fitting of 
many integrations, we find 

no fti 0.91 X lO^cmf ) (— o) ( ^) (^) . (16) 

VO.iy Vl030Gcm^/ \ M ) \m) ^ ' 

Figure 7 shows as an example the dependence of rto on ^. The turn-over radius is less 
than but in all cases close to r^i- Thus, equation (16) is compatible with equation (15). 
Our radius rto has a role similar to that of the Alfven radius rA of Ghosh and Lamb 
(1978, 1979a,b). The dependences we find of rto on M, and M are close to those of 
rA- However, our analysis shows an important dependence on aD which is proportional 
to the magnetic diffusivity of the disk. Furthermore, our rto is smaller than rA by a factor 
of about {aD/12f-^, which is ^ 0.24 for aD = 0.1. 

The turn-over radius is important in the respect that the inward angular momentum 
fiux carried by the disk (MF) is Mouto^to because dw/dr = at r = rto- [Note that we 

14 



do not consider here the possibihty of significant poloidal current flow along the open field 
lines extending from the polar caps of the star (S^ > for 2; > in figure 3) which would 
act to remove angular momentum from the star to infinity.] From numerical integrations, 
we find to a good approximation that ooto = (GM/rf^)^ , so that the infiux of angular 
momentum to the star is MFto = M{GMrto)^ ■ Thus, the rate of increase of the star's 
angular momentum J is: 

^ = M{GMno)--. (17) 

With / the moment of inertia of the star, J = /a;*, and dJ/dt = {dI/dM)M(jj^+I{doj^/dt). 
For the situation of interest, the term proportional to dl/dM is negligible because rto is 
much larger than the star's radius. Thus, we have "spin-up" of the star, 

dw^ M{GMrto)^ 

-dt = 1 ' ^''"^ 

or 

-— ~ 5.8x10 -[Ys)[t^J[ I JKm^J [w^J ' ^^^^^ 

where P = 27r/a;* is the pulsar period. Because the total accretion luminosity L = 
GMM/r^ with r* the star's radius, P ex _p2^o.85^o.285 f^^. constant aD,r^, and M. 

Our outflow solutions are consistent only in the case where Uto = {GM/r^^)'^ > cu*. 
In this case, with r decreasing from rto, the rotation rate of the disk uj{r) decreases and 
approaches the rotation rate of the magnetosphere (~ a;*) from above where it matches onto 
the magnetospheric rotation through a radially thin turbulent boundary layer. We find no 
consistent stationary solutions when uj{r) decreases below a;*. The condition uto > ^* is 
the same as 

.i.5xl0«cm(— j (-j , (19) 

where Vcr is commonly referred to as the "co-rotation radius". Note that Vcr decreases 
during spin-up. Thus, 

Tcr approaches r^o if rto does not decrease rapidly. 
F. Magnetic Braking of Star by Disk; Spin-Up or Spin-Down {vto > ^cr) 
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For rto > fcri the B field in the outer part of the disk (r » Vcr) is open as discussed in 
§IIA, but there are no consistent stationary MHD outflows, {B^)h = 0, (3r = hu/ijt -C 1, 
and the Ohmic dissipation is small compared with that due to viscosity. Thus, the outer 
part of the disk obeys essentially the equations of Shakura and Sunyacv (1973) for a 
Keplerian disk with F = const. = Foo in equation (6) undetermined. We show below 
that Foo can be determined self-consistently by considering the region near where the 
Keplerian accretion flow is brought into co-rotation with the star. Because there are no 
outflows, the angular momentum influx to the star is MF^ which is equal to dJ^/dt. 

In the region of the disk where u{r) is close to w*, that is, where r ~ rcr, \{Bcp)h/Bz\ 
given by equation (2) is not much larger than unity. In this region, closed but twisted held 
lines link the star and the disk. The twist of the fleld acts to remove (or add) angular 
momentum from (to) the disk if B^B^ > (or < 0). The angular momentum removed 
from the disk is deposited on the star via the B fleld. In this region of the disk the key 
equations are (2) and (6) which we rewrite as 

d(jj U , ^ n, , 

^ = (20a) 

^.^(,_,^)(!^^)^(^). (206) 
dr Dvt V M / \rmax^ 

Here, H{x) is a Heaviside function such that H{x) = 1 for \x\ < 1 and H{x) = for 

kl > 1; 

is a measure of the fleld twist; Tmax = const, is the maximum value of the twist; and 

F(r) = r^oj + — . (20d) 

u dr 

We consider Tmax to be a universal constant. For |t| > Tmax the fleld conflguration becomes 
open. In equations (18), Bz is assumed to be the star's dipole fleld. The other equations of 
§IIB arc still needed. Note that both the viscous and Ohmic heating must be retained in 
the energy equation (7). The viscous dissipation is dominant for r > r^r, while the Ohmic 
dissipation dominates for r < rcr- Equations (5) and (8b) can be combined to give u as 
a function of r, T, and u. In turn, this expression for u can be combined with the energy 
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conservation equation (7) to give both u and T as functions of r,uj, and duj/dr if viscous 
heating dominates, or as functions of r and uj if Ohmic heating dominates. Furthermore, 
Vt = {2/3)aCsh in equation (20) can be derived from T and h/r from equation (9) using 
Pr from equation (8b), /3<^ = — r, and e from and u. 

We have solved the equations (5, 7-9, & 20) by numerical integration starting from 
the outside, at a distance where r > Tmax and F = const. = Fqo, and integrating inward 
through Tcr- The fact that the value of is not known a priori points to the use of 
a "shooting method" for its determination. Using this approach, we find that there is 
in general a unique value of Fqo, denoted F^, such that the solution for u>{r) smoothly 
approaches a;* for r decreasing below Vcr- If -P'oo is smaller than F^, then a;(r) follows the 
Keplerian law as r decreases below Vcr- On the other hand, if F^o is larger than F^, then 
u){r) goes through a maximum and decreases rapidly to values much less than co^. The 
only physical solution is that with F^q — F^. Hereafter, the zero superscript on F^o is 
implicit. 

Figure 8 shows the radial dependences of the main physical quantities for a case of 
magnetic braking. The behavior of equations (20) can be understood qualitatively by 
noting that if uj{r) were to decrease linearly through rcr, then the radial width of the 
region of braking, dF/dr < where u! < lv^, would equal that the width of the region of 
acceleration, dF/dr > where a; > a;*. Consequently, F would change by only a small 
fractional amount. However, if within the braking region F{r) attains a value equal to 
a;r^, then duj/dr = at this point and duj/dr remains small and a; ~ a;* for smaller r. 
A rough estimation gives the radial width of the braking region as Ar ~ TmaxF)Cs/ujK, 
where Cg is based on the temperature just outside of this region, and the jump in F as 
AF ~ wAr(r2s2/M). 

Figure 9 shows the dependence of Fqo on ji and on M for sequences of cases where 
Tcr is assumed comparable to rto- This is clearly not a necessary choice; that is, rto could 
differ significantly from Tcr- However, with = "Tto, Figure 9 allows a direct comparison 
of the angular momentum influx to the star (MFx) for the case of outflows (rto ~ ^cr) 
with the influx (MFoo) in the case of no outflows and magnetic braking {rto > i^cr)- Figure 
9 suggests that a bimodal behavior can occur where the star switches between spin-up and 

17 



spin-down, owing to some extraneous perturbation. 
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III. DISCUSSION 

This work has made a new investigation of the problem of disk accretion of matter 
onto a rotating star with an ahgned dipole magnetic field. This is a significant idealization 
in that the relevant objects are mis-aligned X-ray pulsars in binary systems. Nevertheless, 
the aligned case is important for understanding the theory and as guide to theory of the 
mis- aligned rotator. 

The important new ideas introduced in this work are: (1) When the angular velocity 
of the star and disk differ substantially, the B field between the star and disk rapidly 
"inflates" (see Figure 5) and this gives rise to regions of open field lines extending from the 
polar caps of the star to infinity and from the disk to infinity as shown in Figure 3. (As in 
earlier models, there is a closed donut shaped magnetosphere surrounding the star where 
the matter rotates at about the angular rate of the star.) (2) The open field line region of 
the disk may have magnetically driven outfiows or winds from its ±2; surfaces, and these 
outflows can have an important back reaction on the disk structure. We flnd that outflows 
occur when the "turn over radius" rto of equation (16) is less than the "co-rotation radius" 
Tcr = {GM/(jj1)^, where a;* is the rotation rate of the star. The solutions we find with 
outflows from the disk give an influx of angular momentum to the star, that is, the star 
spins up. Thus, a system initially with rto < Vcr evolves with Vcr decreasing towards rto 
if rto does not decrease rapidly. A necessary condition for MHD outflows is that the slow 
magnetosonic point of the outflow be located relatively close to the disk surface, which is 
possible if the angular rotation rate of the disk uj is less than the Keplerian rate by a small 
fractional amount and f3r = {Br)h/Bz = 0(1). If u; is signifcantly less than the Keplerian 
value, MHD outflows do not occur. Under these conditions the angular momentum flux in 
the disk is carried mainly by the bulk transport of the accreting matter and it is inward; 
that is, the viscous transport is small. 

In our treatment of the disk, we assume the "a" turbulent viscosity model of Shakura 
and Sunyaev (1973), and, further, that the magnetic diffusivity of the disk is of the order 
of the turbulent viscosity (Bisnovatyi-Kogan and Ruzmaikin 1976; Parker 1979). This 
magnetic diffusivity is smaller than the anomalous diffusivity invoked by Ghosh and Lamb 
(1978, 1979a, b) by a factor of order 10^. Also, in contrast with the model of Ghosh and 
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Lamb, accreting matter inside of r^o or Vcr (if fcr < Tfo) remains in a disk. Funnelling of 
matter onto the star's magnetic poles may result from from mis-alignment of the stellar 
dipole field. 

When rto > Tcr, we argue that outflows do not occur but rather that there is "magnetic 
braking" of the accretion flow in the vicinity of Vcr- The braking solutions we flnd can 
give spin-up or spin-down (or no spin change) of the star depending mainly on the star's 
magnetic moment and the mass accretion rate. In contrast with the case of outflows, 
the angular momentum flux in the disk is due to both the bulk transport (inward) and 
the viscous transport (outward). For braking solutions with rcr comparable with rto, we 
flnd spin-down of the star for large values of the star's magnetic moment and/or large 
values of the mass accretion rate. Thus, for a system with rcr comparable to rto, bimodal 
behavior is possible where the system flips between spin-up (rto < Tcr) ^ind spin-down. 
This transistion could be triggered by a number of factors, for example, intermittency of 
a, -B-fleld flux introduced from the companion star, or variation in the mass accretion rate, 
owing to the dependence of rto on a, fi, and M. The pulsar GX 1 + 4 shown in Figure 
1 illustrates the bimodal behavior observed for some objects. Future work is planned to 
address the time-dependent accretion onto an aligned dipole rotator including magnetically 
driven outflows from the disk. 
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FIGURE CAPTIONS 



Figure 1. The figure shows the pulse period history for pulsars Her X-1 (a) and GX 1+4 
(b). The first shows mainly spin-up behavior with P = —3.35 x 10~^s/yr. Characteristic 
"wavy" fluctuations of pulse period, with episodes of spin-down, P = 0, are alsoobserved. 
The second pulsar shows steady spin-up with P = — 2.55s/yr (GINGA observations, Bild- 
sten 1993) and then a change to spin down, P = -|-1.5s/yr, and later P = -|-2.51s/yr 
(BATSE observations, Bildsten 1993). 

Figure 2. The flgure shows the model of Ghosh and Lamb (1978). The existence of a 
broad transition zone with closed magnetic field lines was suggested in this paper. 

Figure 3. A schematic drawing of the magnetic field configuration considered in this work. 
The magnetosphere consists of an inner part, where the magnetic field lines are closed and 
an outer part where the field lines are open. 

Figure 4- The figure shows two poloidal field lines used in the derivation of equation (2). 

Figure 5. The figure shows perspective views of illustrative field-line loops in the force-free 
corona of a Keplerian disk. In the bottom panel the difference in the angular displacement 
of the foot points A0 is zero, while the displacement increases going to the middle panel 
(Acf) 80°), and then to the top panel (Acj) 320°). These field lines were determined 
using the equations of Newman, Newman and Lovelace (1992). 

Figure 6. The figure shows the radial dependences of the main physical quantities for a 
case where there are a MHD outflows from the ±2; surfaces of the disk and where the 
pulsar spins up. In the plot of the disk's rotation rate, rto is the "turn over radius" where 
du>/dr — 0. In the plot of F, the radial scale is different from the other plots. In this 
plot rwi is the inner radius of the MHD outflow. We have taken q; = 0.1, /5<^ = 1,D = 1, 
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jj, = 10^°Gcm^, M = 10^'^g/s, and M = Mq. The values of co and ujk are in units of 
(GM/rl)^, Fo = {GMto)^, where = 0.91 x lO^cm. The quantity Tgs is the Shakura 
and Sunyaev (1973) value. 

Figure 7. The figure shows the dependence of the turn-over radius r^o on jj. for a case 
where a = 0.1, (3^ = 1, D = 1, M = lO^V^, and M = Mq. 

Figure 8. The figure shows the radial dependences of the main physical quantities for a 
case where there is no outfiow, but instead magnetic braking of the accretion fiow. For 
the case shown, the infiux of angular momentum MFqo > 0, and thus the star spins up. 
For this plot, a = 0.1, D = I, r^ax = 5, = lO^^G cm^, M = lO^Vs, M = Mq, 
r^r = 0.91 X lO^cm, and Fq = (GMrcr)^ ■ For these values rcr = rto- The quantity Tgg is 
the Shakura and Sunyaev (1973) value. 

Figure 9. The figure shows the dependence of Fqo and Fk = (GMrcr)^ on the pulsar 
period P assuming rcr = Ho for M = lO^^g/s and 4 x 10^^ g/s. The other parameters 
have been taken to be a = 0.1, D = 1, Tmax = 5, and M = Mq. Thus, each M = const. 
curve corresponds to different magnetic moments as oc P^-^'^ . F^o and Fk are measured 
in units of (GMto)^ with To = lO^cm. This plot allows a comparison of the angular 
momentum infiiix to the star [MFk) in the case of outflows and spin- up {rto < i^cr) with 
the case of magnetic braking where the influx is MF^{rto > frc)- 
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